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Abstract 

A scenario for measuring the predicted processes of vacuum elastic and inelastic photon-photon 
scattering with modern lasers is investigated. Numbers of measurable scattered photons are cal- 
culated for the collision of two, Gaussian- focused, pulsed lasers. We show that a single 10 PW 
optical laser beam split into two counter-propagating pulses is sufficient for measuring the elastic 
process. Moreover, when these pulses are sub-cycle, our results suggest the inelastic process should 
be measurable too. 
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I. INTRODUCTION 



Quantum electrodynamics is commonly regarded to be a fantastically successful theory 



whose accuracy has been tested to one part in 10 for free electrons 
in 10^ for bound electrons 



and one part 



2[. However, among its several predictions that have yet to be 
confirmed, is the nature of electromagnetic interaction with the quantised vacuum. Already 
with the pioneering work of Sauter, Heisenberg and Euler js], 0], it was clear that quantum 
mechanics predicts how particles traversing the classically empty space of the vacuum can 
interfere with ephemeral "virtual" quantum states, whose lifetimes are of durations permit- 
ted by the uncertainty relation. Virtual electron-positron pairs, can in principle, be po- 
larised by an external electromagnetic field, thus introducing non-linearities into Maxwell's 
equations, which break the familiar principle of superposition of electromagnetic waves in 
vacuum. Photons from multiple, vacuum-polarising sources, can then become coupled on 
the common point of interaction of the polarised virtual pairs. This process is predicted to 



manifest itse 
one another 



f in a variety of ways such as in a phase shift in intense laser beams crossing 
sl, in a frequency shift of a photon propagating in an intense laser 61, in po- 



larisation effects in crossing lasers such as vacuum birefringence and dichroism 







9|, where 



ideas have already found an applied formulation 10|, in dispersion effects such as vacuum 



diffraction 
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12| and also in vacuum high harmonic generation 13|. The typical scale for 



such "refractive" vacuum polarisation effects, where no pair-creation takes place, is given 
by the critical field strength required to ionise a virtual electron-positron pair, namely the 
pair-creation scale of -E'er = m'^c^/eh = 1.3 x 10^^ Vcm~^ or an equivalent critical intensity 
of Icr = 2.3 X 10^^ Wcm~^, where m and — e < are the mass and charge of an electron 
respectively. Although this intensity lies some seven orders of magnitude above the record 
high produced by a laser [14], recent progress at facilities such as the ongoing 10 PW up- 
grade to the Vulcan laser 15| as well as proposals for next generation lasers HiPER and ELI 
aiming at three to four orders of magnitude less than critical, will put the experimental ver- 
ification of these long-predicted non-linear vacuum polarisation effects finally within reach. 
This therefore motivates more realistic quantitative predictions. 



In the current paper, we focus on the phenomenon of photon-photon scattering, which 
can either be elastic in the sense of a diffractive effect, or inelastic, in the sense of four- 
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wave mixing, allowing the frequency of one field to be shifted up or down in multiples 
of the frequency of the others. When all external fields have the same frequency, four- 
wave mixing is then equivalent to lowest order vacuum high-harmonic generation. As an 
elastic process, numbers of scattered photons have been calculated in the passage of one 
monochromatic Gaussian laser beam through another |l6|], as well as in so-called single- and 
"double-slit" set-ups 0, [lH, where a probe Gaussian beam meets two other intense ones. 
Inelastic photon-photon scattering has been investigated theoreticall y as a four-wave mixing 
process using TEio and TEqi modes in a superconducting cavity [l7|, in the collision of 
three, perpendicular, plane-waves U and as generating odd harmonics involving a single, 
spatially- focused monochromatic wave [19]. By incorporating both the pulsed and spatially- 
focused nature of modern high-intensity laser beams, we perform a more accurate calculation 
of the signal of the elastic scattering process. We thereby investigate the robustness of 
the effect with a more detailed calculation than hitherto performed, including dependency 
on beam collision angle, impact parameter (lateral beam separation), longitudinal phase 
difference (through lag) and pulse duration (finite beam length). Inclusion of four- wave 
mixing terms with a pulsed set-up allows us, moreover, to determine the possibility of 
measuring inelastic photon-photon scattering when a single 10 PW beam is split into two 
counter-propagating sub-cycle pulses. In what follows, we work in Gaussian cgs units (fine- 
structure constant a = e^), with h = c = ^tteq = 1, unless explicit units denote otherwise. 



II. SCENARIO CONSIDERED 

In order to analyse the collision of two laser pulses, several collision parameters have been 
included. The envisaged scenario is shown in Fig. [H in addition to which, lateral and 
temporal centring and carrier envelope phase appear in the analytical set-up. Spatial focus- 
ing and temporal pulse shape are present in taking the leading order spatial and temporal 



terms of the Gaussian beam solution to Maxwell's equations (see e.g. 20]). These approxi- 
mations neglect terms of the order 0{wcfl/yr,c) and 0{1/ijJcTc) respectively, where c G {a,b} 
is used throughout for beams a and b, the minimum beam waist is Wcfi, Rayleigh length 
yr,c = ^cw1q/2, beam frequency Uc and full-width-half-max pulse duration tfwhm related to 



r via rv21n2 = tfwhm- The condition uJcTc ^ 1 limits the minimum pulse duration that 
can be consistently considered in our analysis. For the electric fields of the two beams £„, 



Eh, we then have: 



Ea{x,y,z,t) =e 
Eb{x,y,z,t) =ih 



Vl + (y'/2/r,a)^ 



: sin [V^, + Ua{t -At + y') - r/„(y')] fa{t - At + y') (1) 



^c(z/) = tan M 



yr,c 



sin [V^b + - + 776(l/)] /6(t - y) 



2 1/2 + ^: 
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(2) 
(3) 



where the co-ordinates (x, y', z') are the same as {x, y, z) rotated anti-clockwise around the x 
axis by an angle 9, with the polarisation e'^ being similarly rotated so that k^- = k'^' ~ 
and led = l^cl = where kc is the beam wavevector, /c describes the pulse shape with 
fc{x) = 6"'^^/'^'=^^ being used, Wc is the beam waist w"^ = w^q(1 + (?//?/r,c)^) dependent on 
transverse co-ordinate, ipc is a constant phase. At is the lag and Ec^ is the field amplitude, 
which satisfies J dt dx dz \E{x,y = 0, z, t)p/(47r) = S, with total beam energy S, or Ec^ = 
l^fTP^^jwci^, for peak beam power Pc,0) where we have already assumed that corrections to 
transversality k A E = B can be neglected, being as they are, of the same order as neglected 
higher-order terms in the spatial Gaussian beam solution to Maxwell's equations. 
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Fig. 1. The envisaged experimental set-up. t^}, refer to the pulse durations in the Gaussian beam 
envelopes e"^*^^-* l'^a,b ^ the beam is displaced from the y-axis by co-ordinates xq, zq^ both beams 
have in general a carrier-envelope phase and the Eq beam lags behind E;, by At. The E;, field is 
incident on the detector. 



We focus on the phenomenon of diffraction and specifically the detection of photons whose 
wave-vectors differ significantly, either in orientation or in magnitude from those of the 
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background lasers. As such, we envisage an array of photosensitive detectors being placed 
some distance away from the collision, y^, along the positive y-axis, much larger than the 
interaction volume (the subscript d refers to quantities on the detector). is then incident 
on this detector. 



III. DERIVATION OF SCATTERED FIELD 

When external electromagnetic fields polarising the vacuum have equivalent photon energies 
much less than the electron mass {hw <^ mc^), their evolution can be well-approximated 
by an effective description in which the vacuum fermion dynamics has been integrated out 



and only photon degrees of freedom remain. The Euler-Heisenberg Lagrangian {4] is an 
effective Lagrangian which includes such fermion dynamics to one-loop order. When the 
field strength is much less than critical {E ^ E^r), the Euler-Heisenberg Lagrangian can be 
well-approximated by its weak-field expansion, which, neglecting derivative terms, is: 

(4) 

The weak-field expansion Eq. (jl]) is depicted in Fig. H] and can be understood as coupling 



O 





+ 

2n 



Fig. 2. The Euler-Heisenberg effective Lagrangian is an integration over the high-energy (fermion) 
degrees of freedom. The external field can be generated by multiple sources, indicated by photons 
in the diagram being with and without crosses. 

the flux of electromagnetic fields from different sources with one another. Extremising Eq. 
dl]) with respect to the photon gauge field returns the wave equation for E and B fields 
modified by the one-loop, weak-field, vacuum current, Jvac: 

V^E - d^E = 47rJ.ac (5) 
Jvac = fv A dtM - V(V ■ P) + dM , (6) 
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where P = 3-E M = — + — B and 



a2 



P= 180^ + 
M = -^^^ [2(E^ - i?^)B - 7(E . B)E] . (8) 

Using the beam transversahty, A = B^ P and M can be written entirely in terms 
of the electric or magnetic field. One can then write the vacuum polarisation as a series 
Pi = XifEj + xlfkiEjEkEi + . . . xfj^-Vi^Ej ■ ■■Em + . . ., where electric susceptibilities x only 
occur at odd orders due to charge-conjugation symmetry (Furry's theorem). Therefore, 
four-, six-, eight-, etc. wave mixing can in principle occur, although each extra order will be 
suppressed by a factor a{E / Ecr)"^ . 

An iterative approach can be used to solve Eq. ([5]), which, since J*^*^-* oc [a(-E'''°V-Ecr)^]E*-°^ 
and a{E / EcrY <^ 1, can be understood as perturbative: 

/t(")i'F'(") "pCO)"!! 
d'x -" >^'='-\ (9) 

where E*^"^ is the n-th order perturbative solution of Eq. ([5]), E^*^) is the zero- field vacuum 
solution, obeyed by the Gaussian beams in vacuum, approximated by E^ and Ej,, J*^") is the 
ra-th iteration of the current occurring on the right hand side of the wave equation, (x^, td) 
are the co-ordinates in the detector plane and tret = td — — x| is the retarded time. By 
making the approximation that: 

Ej^)(x„t,) = / d'x ' ^'*=*- « EC'), (10) 

J |x-Xd| 

for all n, the resultant electric field can be well approximated by E = E*^^) = E*^°) + E^, 
Erf = E^°\ i.e. the zero-field vacuum solution plus the lowest order "diffracted field." 

By substituting E*^*^) = Eq + Ef, in Eq. ([6]) and Eq. ([9]), and by enforcing the assumption 



that the dimensions of the interaction vo 
co-ordinates, following similar steps to 

2 . 12 



jime are much smaller than the typical detector 
8|, mi, one arrives at: 



i=i ^ 

^) = -^^J^}l E ^^EahEl^iu)^^^!.,, (12) 



where It,j and Iu)j are integrals over the interaction volume, given in Eqs. (1A.3P and (lA.4p . 
Td = |X(i| is the detector distance, with Aj G {—2, —1, 1, 2}, Bj, Tj G {1,2} coefficients given 
in the appendix in Tab. HI r,- = [Bj/t^ + Tj/r^]''^/'^, uj = l3jUa + IjOJb-, i>j = l3ji>a + Iji'b, 
(3j, 7j G {—2, —1, 0, 1, 2} also given in Tab. Hand /(j) = 1 for j <= 6, otherwise = 2 and 
Vi 2 are the diffracted field polarisation vectors given in Eqs. (lA.ip and (IA.2p . 



Splitting the plane- wave part of the input fields E^, Ef, into positive and negative frequen- 
cies, the twelve terms in Eqs. (fTT]) and ( !T2|) are produced, corresponding to the six possible 
orientations of the currents connected by the effective vertex in Fig. |21 As the interaction 
contains terms of the order 0[(Ea + E^)'^], and as the purely cubic terms E^f^ necessar- 
ily disappear (both electromagnetic invariants are zero for the individual Gaussian beams, 
transverse in this approximation), for an incident current of frequency Ub, the resultant sig- 
nal can have a frequency Ua, cOb, ^b ± 2wa, 2ujb ± w^, corresponding to the two beams' elastic 
and inelastic components respectively. The diffracted field polarisation vectors vi^2 appear 
as geometrical factors and from their definition in Eqs. (lA.ip and (IA.2p . one can see that 
on the detector {yd > 0), the Ua and Ua ± 2ijJb signal from pulse a are strongly suppressed, 
as would be expected as the E^ pulse travels from the interaction region away from the de- 
tector. After a further analytical integration in x, the remaining two-dimensional integrals 
from Eqs. (IA.3P and (IA.4P were then evaluated numerically in C++, partly using the GSL 
library 21 1. 



One can interpret the classical field incident on the detector as being composed of a to- 
tal number of photons Nt by dividing its total energy by the photon energy so that Nt = 
J^^dudxddzdlt{uj,rd)/\uj\, where the total spectral density It{uj,rd) = \Et{^,rd)\'^ /Stc'^ = 
\'Eb{u,rd) + E(i(a;, r^) p/Svr^ (Ea{oJ,rd) is taken to be zero in the current beam set-up) and 
where yd is taken large enough that the surface perpendicular to the Poynting vector can be 
well approximated as being flat. Although the spectral density extends to negative frequen- 
cies, it is consistent to interpret the differential number of photons as this divided by the 
absolute frequency because the total energy is the integration over all frequencies and all 
energy is carried by positive-frequency photons (see also I22I] on this point). We then calcu- 
late the number of "accessible" photons that fall on the detector plane, by integrating over 
the annulus that satisfies dNt{xd,Zd) — 100 dNb{xd, Zd) > for dNi = J^^du Ii{ijj,rd) /\uj\, 
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i G {a, b, d, t}. 



IV. ELASTIC PHOTON-PHOTON SCATTERING 

Current and next generation high intensity lasers will typically produce pulses with many 
optical cycles and so unless some resonance condition is fulfilled, one would expect the elastic 
cross-section, where incident and outgoing spectra have the same form, to be the largest. 
By "elastic," we are therefore referring to terms in with equal incoming and outgoing fre- 
quencies. As an analytical test of our expressions, we can reproduce the electric field derived 
for the three-beam, double-slit cases given in j3, llH when the separation of the slits is sent 



to zero - the two-beam limit (this limit was calculated for ll| in 23|), which we label E^, E^ 



d 



referring to head-on and perpendicular collisions respectively. By taking xq = zq = = 
and the limit r^^b — cxd in Eq. ( ITT]) , with ^ = 0, we recover E^ as given in 23|], and with 



6 = 7t/2, we recover E^ as given in p]. As a numerical test of our expressions, we can com- 
pare Idit, Tfi) = |Erf(t, Trf) ■ Edit, rrf)p/47r in the case /^(O, r^) with results using the single-slit 
version of E|j. The equivalent parameters are = 0.8 /im, Aj, = 0.527 um, w^.o = 0.8 /xm, 
Wb,o = 290 /im. Pa = 50 PW, Pb = 20 TW, as the field strengths in 0, Q were calculated 
using a conservative form of the beam intensity with power per unit area for an area vtw^Q) 
rather than the 7rw^Q/2 which is manifest from an integration of the intensity of a Gaussian 
beam over the transverse plane. In order to obtain the agreement between E^^ and E^ shown 
in Fig. [31 Ta,b had to be set to around 10"^ fs, which is unexpectedly large compared to the 
pulse durations considered in those references (r^ = 30 fs, = 100 fs). We will elaborate the 
non-trivial dependency of Id on pulse duration, which explains why most of the difference 
between Id and disappears already at r = 10^ fs. When the number of accessible photons 
was calculated for the pulsed system with 6 = 0.1 and the same durations as suggested in 



231], the number of photons A^^ also fell from the estimated value of around 36 to around 0.4. 



To illuminate the two orders of magnitude difference in A^^ for these parameters, the inte- 
grand for Ed was reduced to the most significant terms for a head-on, elastic collision and 
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Fig. 3. Numerical comparison of Id with results in [23] (denoted I^^)- Plotted is log^Q of the 
absolute relative difference in for, from top to bottom, Ta = ti, = 100 fs, 1000 fs and 20000 fs. 
Between the dotted lines, the background dominates (/;, » Id)- 



evaluated independently in Mathematica. The simplified expression Nd^r) was then 



dpdPd 

max 

2 / dpdPd 



°° du 
-oo 27r 
r 



dy 



V{y,u,Ub,T, pd/vd) 



dy lim lim V{y,u,Ub,T, pd/vd) 



?>m-nydlcri 



2 ^^6,0 1 + {y/yr,h) 



(13) 
(14) 



exp 



<o l+{y/yr,a)■'^ 3 3r2 12 ^ 



4r^[2 + 



""M l + (2//2/r,6)' 



where = {2/wIq + 1/w^q)~"^, = (x^ + z1)/r'\ and A*"^ is equivalent to the expression 
leading from E^. The dependence on pdj^d of these two expressions is shown in Fig. Hl^a), 
where it can be seen that the monochromatic A^^ is much larger and more sharply peaked 
in the centre of the detector. After integrating between the relevant annulus of p~^^^^ jvd = 
0.0032 and p^"""""^ /vd = 0.03, this independent test then gives Ndir) = 0.33 and N^l = 37.0. 
The corresponding values for a circular detector of radius 15 cm from the full expression 
are Ndir) = 0.33 and A"^ = 38.0, supporting the two-orders of magnitude difference and 
the claim that Eqs. fll3H15l) incorporate the main physics. By plotting the exponential 
dependency on Pd/rd, which is integrated over to acquire the expected number of photons 
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Fig. 4. Left-hand plot (a): log^^Q A''^ (upper curve) and log^^Q A''^ (lower curve) plotted against p^. 
Right-hand plot (b): a log-log-plot of Re [l{y)] (solid line) with the first, second and third dashed 
lines corresponding to lim^^ (,_^oo ^y), \l{y)\ and limy,,^^oo /(?/)• 



in the monochromatic case, A*"^: 



l{y) = exp 



^ Pd(i + iy/yr,a) 



(16) 



^V^^ Ub.o^ i+{y/yr.b?i 

we observe the interesting behaviour shown in Fig. SJ^b). We first note that the decay is 
not purely exponential, but has two important length scales: Whfi and which, in the 
limit of being infinitely large, correspond to the first and third dashed curves in Fig. IH^b). 
The inclusion of these extra longitudinal length scales in 23(|, which cannot contribute to 
photon scattering when the finite pulse length is taken into consideration, then explains the 
discrepancy in the values of N^ij) and A*"^. Only the region of the pulses within a distance r 
around their maxima in the longitudinal direction can efficiently contribute to the scattering 
process, with the rest of the pulse being damped by its Gaussian shape. The finite length 
of laser pulses probing vacuum photon-photon scattering can then only be neglected, when 
the duration r is the largest longitudinal length scale. In the limit r ^ yr^b in the full ex- 
pression for Ed in Eqs. (ITT]) and (IT^ . the scaling Ndij) oc r of [23| is recovered, supporting 
this statement (this will also be apparent from Fig. [S]). Furthermore, the results of 23 1 
are expected to remain valid in the case tiwIq/ XbTa^b < 1, so for more focused and longer 
wavelength probe beams as well as for longer pulses. Indeed for the parameters quoted, that 
the effect would be two orders of magnitude weaker is in no way prohibitive to conducting 
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such experiments. For example in |23| . the 
only around Ip ~ 10^^ Wcm~^, but as A^^^ oc Ip 
the probe beam more (if Whfi is set to 60 /im, 
the single-slit and 4.5 in the double-slit case) 
10 TW), to which A^^^ is proportional. 



intensity of the probe beam was taken to be 
, the shortfall could be made up by focusing 
Nd increases approximately by a factor 7 in 
or increasing the power of the probe (from 



The current treatment also allows for the two lasers to be equally strong and we consider 
the more experimentally-accessible situation of having a single laser, split into two colliding 
pulses, both focused to ultra-high intensities. Since scales with E^q if we keep the 
power of the laser constant (-Eco = 2a/2^^/wc,o), for each term, the optimal division of the 
total power between the beams is: 

Pa^O Pj j-f Pj 7") / i~r\ 

= TT'i Pafl = — -jrPtfi- (17) 

For base parameters similar to that of the Vulcan laser |l5] Xa = Xb = 0.91 /im, Ta = = 
30 fs. Pa = 5 PW, Pb = 5 PW, with Wa^ = 0.91 /im, Wbfi = 100 fim, = = x a summary 
of the dependency of on several variables is given in Fig. |5l We will comment on the plots 
sequentially, in which solid lines represent what one could intuitively expect, as explained in 
the following. Starting from the right-hand side of the first plot and moving in the direction 
of falling Wbfl, we see Nd^Wbfi) increases approximately as oc w^q, indicated by the solid 
line. Since Nd for such a set-up is proportional to -E^c ^^"^ since this is inversely propor- 
tional to the area of focusing, the dependency on oc w^^ is as expected. Deviation occurs 
when a maximum is reached (see e.g. for details), beyond which Nd{wbfi) falls rapidly 
as the background from gradually covers the entire detector, leaving no signal. The de- 
pendency on beam-separation Nd{xo) is also intuitive and seen to have excellent agreement 
with a Gaussian, normahsed in height, with a width of Wbfi/2 (exp(— 2a;o/w^o))- Simply 
by integrating the transverse Gaussian distributions of the two beams, and then squaring 
{Nd oc |Edp), one arrives at this dependency. The third plot of Nai^X) {X = Xa = Xb) is 
a log-log plot where the dependency begins for small X as Nd ^ X~^ but then for larger 
values tends to A^^ ^ X~^'^. This is shown by all the points lying between these two solid 
lines. Since the power of each beam is inversely proportional to wavelength, and since the 
Nd oc P^QPbfi, one would expect at least a dependency of Nd{X) ~ A^'^. In contrast, the 
dependency of Nd{T) can be straightforwardly derived. For Ta = u = r, one notes that 
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when r ^ Wbfl, yr,b, the interaction volume in beam propagation direction is governed by 
the Gaussian pulse shape. Further noting that A^^ essentially involves a double integration 
on longitudinal beam co-ordinate (through taking the mod-squared), as well as an integral 
over t, the dependency N^ir) oc appears, which shows excellent agreement for small r 
with the full numerical integration, displayed by the on the log- log plot of Nii{t) in the 
fourth figure. The larger r is for r > Wbfi, the more the decay along the beam propagation 
axis is described by focusing rather than pulse terms. For large enough r. Id depends only 
on focusing terms and since the yield is acquired from an integration over time, we have 
Nd{T) oc r and this transition can be seen by the second, linear, fit line for large r in the 
figure. An estimation of the dependency of Nii{6) oc (1 + cos^)^ on beam intersection angle 
comes from the geometrical factor in vi oc (1 + cos^^), which must be squared and gives the 
approximate agreement shown in the fifth plot. For small angles Nd{6) cc 1 — 6"^ /2, making 
the dependency relatively weak for near head-on collisions (A^^^ remains at 90 % of its value 
up to 6 ^ 7r/7). The final plot of Nd{At) closely resembles a Gaussian with width 9r and 
so for this set-up, A^^ is relatively insensitive to lag. 



One strategy to increase the number of diffracted photons would be to use higher-harmonics 
of the probe laser. If the same parameters as in Fig. [2] are used, for a collision angle of 
6 = 0.1, assuming a 40% reduction in energy due to generating the second harmonic, A"^ ^ 4. 
If this process could be repeated to generate the fourth-harmonic, with a 16% reduction, 
Nd ~ 13. As previously argued in 11[, such numbers of scattered photons should allow 
detection in experiment. A discussion of sources of background noise and why they can be 
effectively neglected is given in 
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23|. 
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Fig. 5. Dependency of the number of measurable elastically diffracted photons on various 
parameters, where parameters held constant take the values Aa = A;, = 0.91 /im, Wafl = 0.91 /xm, 
Wb,0 = 100 /xm, Ta = n = 30 fs, Pa = 5 PW, = 5 PW, Sa = h = X. 
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V. INELASTIC PHOTON-PHOTON SCATTERING (FOUR- WAVE MIXING) 



When considering the possible frequencies of the resultant current, conservation of energy 
and linear momentum leads one to the equations: 

u = sgn(/3j)K,i + Sii3^\2UJa,2] + sgn(7j)[wfe,i + (^17, 12^^,2] (18) 
+ sgn{'yj)[uJb,i cos 6*6,1 + <5|7,|2C^b,2 cos 6*6,2] (19) 

+ sgn(7j)[wb,i sin6'b,i + 5\j^\2UJb,2 sin 6*6,2], (20) 

where cu is the frequency of the resultant current, sgn(a;) returns the sign of x with sgn(O) = 0, 
and 6'{a,6},{i,2} are the angles the currents make with y. Therefore, detection co-ordinate, 
focusing and harmonic order are already linked at this stage. It turns out to be difficult to 
satisfy these conditions simultaneously with just two laser beams and a fixed observation 
angle. For example, if we take /3j = 2, 7^ = 1, with C0a,i = uja,2 = i^a, 1^6,1 = i^b,2 = 
for simplicity and a more-or-less head-on collision of the lasers, so 6'a,{i,2} is approximately 
equal to vr — 6*6,(1,2} and 6'6,{i,2} is small, then, to first order, from Eqs. (fT8|) and (fT9|) we have 
oj = 2uja + oJh and ujyd/rd ~ — 2u;a + ujh- Since ydl^d ~ 1 on the detector, the contribution 
from this term can therefore only be satisfied by a small range of frequencies around = 0, 
which are not typically populated in the spectrum of E^. The energy- momentum conditions 
Eqs. f lT8H20|) can be most easily seen occurring in the exponent of the integral /^j Eq. flA.4p . 
where they appear as frequencies of plane waves to be integrated over in y, 2, becoming 
Gaussian-like after integration. The larger the deviation from these conditions, the higher 
the frequency of oscillation to be integrated over, the more exponentially small the resulting 
amplitude, typical for evanescent waves. 

We investigated the ansatz that for short enough pulses, the bandwidth of the two lasers 
becomes wide enough that Eqs. f lT8H20|) can be fulfilled simultaneously for a measurable 
amount of photons. Essentially, for this four-wave interaction, three different photon energies 
can be supplied by two lasers. To make this statement explicit, instead of using a temporal 
envelope, we can consider building the pulses in the frequency domain: 

/oo 
rfa;eE~(x,t,a;e)^7(a;c,Wc,o), (21) 
-00 

14 



where E™°'^°(x, t, Wc) is the electric field of a monochromatic Gaussian beam, frequency Uc 

and g{(jJc,0Jcfi) is the spectral density of the pulse E^(x, t), with peak frequency w^o- Then 

B - r 

due to our interaction being cubic in the fields {Eia\ E^^), the integration over this current 
in the frequency domain, Eq. ( 1A.4I) . would include three extra integrations over frequency 

/ duadubducg{oJa,uJa,o) g{uJb,u}b,o) g{uJc,6B,2U}a,o + 5r,2Wfe,o)^(^ - ^i)' where the final delta 
function appears explicitly from an integration over t. Here it is apparent that due to the 
finite bandwidth, in general, three different energies enter the effective vertex in Fig. [2] 
from the two lasers. If the spectrum is taken to be Gaussian g{uc,u}cfi) = exp[— — 
|^c,o)^Tc /4]rc/2Y^ we have, setting = without loss of generality: 

2 2 2 

E' x,t =e, / dooc ' , = e +"'<o«^ + c. c. 



oo 

2 , ,2 



■ sm 



/ \ -1 y ^cov X + z 

^c,o{y - t) + tan — 



+ h. o. t.. 



Vr 2 y'^ + y^.- 

where the remaining terms are of the same order as those neglected in the Gaussian beam 
solution. Therefore the use of a Gaussian temporal envelope in E^ and E;, (Eqs. ([1]) and 
([2])), is equivalent to integrating over three different photon frequencies from the external 
fields in the interaction. 



When xo = zo = At = 6 = 0, yr = yr,a = yr,b and pl/rj = {xj + zj)/rj is small, dNd{xd, Zd) 
can be approximated analytically. We can write dNd{xd, Zd) = J2l\=i dN^'^{xd, Zd) and 
demonstrate this analysis by concentrating on a single term A^J"' for convenience (the full 
expression is given in Eq. (]A.9P ). One can show: 



dNf{xd,Zd) ^ 



7r2 



90 El 16 STd 



2 , ,2^2 ^2 ^2 

2~t-^+ 

g 11 



dUj\Uj\ e d ^ ''1' 99 , (22) 

) 

where s = l/{Bj + Tj{wafi/wbfiY), Tgg = Tg/fg, = {Bg/rl - Tg/r'^)-^ and Ugg = Ug/ug, 

g = PqUJa — 'JqUJb, uudcr the condition T'^/y^ ^ 1, for = r^[l — t^^]~^ and where a 



condition on co: |[w(rgg — yd/fd) + ojg — ujgT'^^'^ ly^\ ^ 1 has been approximated by taking 
the upper limit of the integration as oo. To simplify the discussion, let Tq = Tb = t. Then 
we can see from Eq. (!22l) that the spectral density for inelastically scattered photons has 
a different shape to the background, namely with a minimum at a; = and two maxima. 
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whose positions for the case Xd = Zd = are u"^ = {'ygUb/2){l ± [1 + 12/(7qa;br)^]-^/^). 
Using a spectral filter, and short enough pulses, this could in principle be used to separate 
the different inelastic scattering signals from each other and the elastically scattered and 
background photons on the detector. Setting Pd = for brevity, the final integral can be 
approximated by: 




1 2 



W 



a,0 



180 El 16 su 



Erf 



V2 



(23) 



It should be noted that v/(j) is identically zero for j > 6 at = yd, and so the frequen- 
cies Ua, 2ijjb ± uja are suppressed, as already argued. The numerical integration of the full 
highly-oscillating integrands was performed using the Filon method, which is an approxima- 
tion to the integral / dtf(t) cos{ut) for asymptotically-large u (see e.g. [24]), used with the 
GNU arbitrary-precision C++ library [25]. Agreement between numerics and analytics for 
'u^afi = U!b,o = 10 /xm, = 1 m, /3g = 2, 7q = 1 is then shown in Fig. |H1 in part corroborating 
our numerical approach. 

_ 0.016 



1 

OH 



0.01 



0.005 



1.5 
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T [fsl 
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Fig. 6. Agreement of the analytical approximation Eq. (I23p with the corresponding numerical 
solution. 



The pulse duration of each laser plays an important role in four-wave mixing. By choosing a 
temporal profile for the beam that is Gaussian, we already have implicitly the lower bound 
r ^ l/cj. As pulse duration and longitudinal co-ordinate are linked, a natural upper bound 
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is also formed for our calculation in the assumption that the diffracted field is smaller than 
the vacuum-polarising fields Eq. f lTO|) . Assuming scattered photons arriving at a point on 
the detector are generated in the centre of the beams' intersection, the integration is ex- 
clusively over regions in which the polarising beams are more intense than the diffracted 
field when r ^ 2wfe oZ/d/Pd, giving l/cu ^ r ^ Sw^ oZ/d/Pd- The lower bound limits our abil- 
ity to assess the importance of the inelastic process. We require a large bandwidth Au/u 
for the inelastically-generated photons to be on-shell, but from the bandwidth theorem, 
Au/u ~ l/wr ^ 1 by our limitation on r. As a consequence, with a two-beam set-up, spec- 
trally separating off the inelastic signal would be experimentally challenging, as this signal is 
generated when the bandwidth of the elastic background overlaps these "inelastic" frequen- 
cies. More promising seems to be to observe the change in due to inelastic scattering 
becoming significant as r is reduced. In Fig. [71 we plot this ratio {Nt — Ne)/Ne against r^, 
where N^, is the number of photons scattered due to when only the elastic terms are included 
in Eq. ( fT2|) . The results suggest that for short enough pulse durations, the inelastic process 
can influence the total number of measured photons substantially. In Fig. [7] the proportion 
reaches over 20%, for a minimum pulse duration of r = 1 fs, equivalent to UaTa ~ 2. This 
could already have been anticipated from Ed{^d,td) in Eq- ffTTI) . including, as it does, a 
pre-factor 4 + {ujTjY. In addition, although the pulse durations are short, assuming again 
40% attenuation each time a second-harmonic is generated from the probe, the total number 
of diffracted photons ranges from 1 to 4 (at = 1, 2 respectively). Although the analysis is 
limited by how small can be consistently made, these results lend support to the ansatz 
that two laser beams with a large bandwith, especially in the laser being probed, can be 
used to measure the effect of the inelastic process. 

In order to further support this ansatz and without being limited by a minimum value of 
the pulse duration, we can consider the simplified case of the collision of two plane waves 
modulated by a sech envelope. 



L n J 

These fields satisfy Maxwell's vacuum equations exactly, therefore removing the limitations 
on conceivable pulse lengths brought about by using a perturbative solution. The analysis 




a 



(24) 




(25) 
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Fig. 7. The increasing importance of the inelastic process with increasing bandwidth. Plotted is 
the proportion of the total number of diffracted photons that are due to inelastic scattering, against 
Ta, for Pb = 10/3 PW, Pa + Pb = 10 PW, \a = 0.91 /um, \b = 0.2275 /im, n = 2 fs, Wa,o = 0.91 /im, 
Wbfl = 50 /im, Sa = (1,0, 0), Eb = (0, 0, 1), 1pa = tpb = 0. 



proceeds just as for the Gaussian case but with the difference that now the fields are not 
bound in the transverse plane. Therefore, in order to avoid a divergence, we only consider 
the resulting P and M to be non-zero up to a finite transverse radius po- It can be shown 
that this curtailing of the interaction region then allows us to integrate over the current Eq. 
dH]) as usual. The diffracted field E^'^'^'^(xrf, w) then becomes: 

(26) 



Abn^m'^rd 8i ^-^ 



1 ^ 



where are geometrical factors as in the Gaussian case Eqs. (lA.ip and (lA.2p . I^'^^ are 
integrals given in Eq. f lA.lOp . the sum over j corresponds to the two terms El^Ebfi and 
EafiEl^ respectively and = has been set for simplicity. Unlike for Gaussian beams, 
the elastic scattering terms cannot be isolated so easily. In order to exemplify the effect of 
the inelastic process however, one can observe how the behaviour of N^^'^^^^ changes as oOaTa 
is reduced to below unity. Deviation from "elastic" behaviour, indicates the importance of 
inelastic scattering. 



The first plot in Fig. [8] depicts the dependence of N^'^'^ on and we notice that for UaTa ^ 2 
{Ta < 0.8 fs), there is indeed a deviation in the behaviour of N^^'^^\ We can take data from 
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a more uniform region Tq > 1 fs and acquire a best-fit polynomial with the boundary 
condition Nji'^^^\Ta = 0) = 0. It turns out that a cubic polynomial fits the calculated 
points well (similar to the Gaussian beam case where Nii{T) oc r^). When the fit parameters 
were calculated for 1 fs < < 2 fs, the goodness-of-fit was tested with a Pearson's chi- 
squared test over 1 fs < < 3 fs and found to support the hypothesis of agreement with a 
probability of over 0.995. When the relative difference of this "elastic" curve from the total 
was calculated, the second plot in Fig. [8] was generated. This clearly demonstrates the new 
behaviour occurring for short pulse durations or equivalently large bandwidths and so further 
supports our initial ansatz that just one beam split into two counter-propagating sub-cycle 
pulses is sufficient for accessing the process of vacuum inelastic photon-photon scattering. 
A suggestion for further work would be to investigate the role of the carrier-envelope phase 
as well chirped frequency. 



0.6 r 




0.001 0.01 0.1 1.0 2.5 ■ 0.1 0.5 1 1.5 2 2.5 



Ta [fs] Ta [fs] 

Fig. 8. On the left, for N^^'^^ira), with Tf, = 1.4 fs, at Xd = 0.1 r^, = 0, = 1 m and 
Aa = = 0.8 fim, the dominant term E'^Ef, has been plotted and the coefficient of the integral 
ignored. On the right is plotted the relative difference between the "elastic" and full behaviour of 



VI. SUMMARY 

In calculating numbers of photons scattered in the collision of two laser beams, we had three 
aims: i) to consider a more realistic set-up of the colliding beams (including a temporal 
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pulse shape, collision angle, lag and lateral separation), which would produce more accurate 
qualitative and quantitative predictions for experiment, ii) to investigate the possibility of 
using a single laser, split into two beams to measure elastic photon-photon scattering and 
iii) to evaluate the ansatz that just two lasers, with sufficiently short pulse durations, can be 
used to measure the process of inelastic photon-photon scattering. The first of these aims has 
been met in Fig. [5] where the dependency on various collision parameters was calculated and 
found consistent with physical reasoning. This led to the second aim, where the inclusion 

ed to two orders of magnitude difference over previous 



of a pulse form and collision angle 
elastic photon scattering estimates [23| (the single-slit limit of Uj). In this more complete 
description, it was shown that when a 10 PW, A = 0.91 nm beam is separated into two 
30 fs Gaussian pulses, incident at an angle 0.1, one could expect approximately 0.7, 4 or 
13 photons, corresponding to the fundamental, second and fourth harmonic of the probe 
respectively (with an assumed loss of 40% per frequency doubling), to be diffracted into 



detectable regions. As argued in [ll|, this could be sufficient for measuring elastic photon- 
photon scattering, here shown using a single 10 PW source. The final aim was partially 
met, first by considering Gaussian pulses, where it was shown that for uoTa < 4 for the more 
intense beam a, the inelastic scattering process increased and became as large as around 20% 
that of the elastic count for uoTa ~ 2. However, for these results to be consistent, UTa ^ 1, so 
the head-on collision of two sech pulses was analysed, for which no such bound applies, where 
it was shown that again, in this different field background, for uTa ~ 2, inelastic scattering 
became important - as large as around 50% that of the elastic one, lending supporting to 
our original ansatz. 
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Appendix A: Integration formulae 
1. Gaussian diffracted field formulae 

Sum coefficients: 
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TABLE I. Sum coefficients that occur in the expressions for and E^, Eqs. (jlip and (I12p 



Diffracted field polarisation vectors: 



Vl 



4 (el ■£,(! + cos 0) - y' . y . e'^) [e'^ (l - (" r r)) + (y' A e'J A r 



7(< -yAet-ef y'Ae) - (y' Ael)( 1 - f- r 



V2 = 4(£l -£,,(1 + cose) -y'-ib y-e'A eJl - r r) ) - {y' Mb) A r 



+?(< ■ y A £6 - • y' A el) [(y A (l - ("^ r r)) + fff, A r 



(A.l) 



(A.2) 



Integration terms: 



,2., , , „ll2 



Z!2Tf ^ + s'^',^2!!,,2 > +^ ^C!2!r„2"°? +i{^aPjy' -^Ij^^bV-Pj tan i(y/y,,a)+7, tan i{j//j/r,6)) 



2{!/^+!/;;_i,) 2(y'2+„2^) 2(i,2+y<!^) 



e 



_2t(£ll4^-£^)+2(r,--^)(4 
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(A.3) 



^^{x^^'z''^) |[(a;-a;o)^ + (^--S())^]+t(/3jC^ay-7j'^6?/)-2/3j tan ^(yV?/r,a)+^7j tan ^{y/yr,h) 
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2. Analytical approximation to dN^ 
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3. Sech diffracted field formulae 



Integration terms: 
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PO /•27r 

dp p d(p 
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